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We realize indirect partial measurement of a transmon qubit in circuit quantum electrodynamics 
by interaction with an ancilla qubit and projective ancilla measurement with a dedicated read- 
out resonator. Accurate control of the interaction and ancilla measurement basis allows tailoring 

the measurement strength and operator. The tradeoff between measurement strength and qubit 
back-action is characterized through the distortion of a qubit Rabi oscillation imposed by ancilla 
measurement in different bases. Combining partial and projective qubit measurements, we provide 
the solid-state demonstration of the correspondence between a non-classical weak value and the 
violation of a Leggett-Garg inequality. 



PACS numbers: 03.67.Lx, 42.50.Dv, 42.50.Pq, 85.25.-j 

Quantum measurement involves a fundamental trade- 
off between information gain and disturbance of the 
measured system that is traceable to uncertainty rela- 
tions [1] . The back-action, or kick-back, is a non-unitary 
process that depends on the measurement result and 
pre-measurement system state. Thought experiments in 
the 1980's unveiled paradoxes where the back-action of 
multiple measurements of one system puts quantum me- 
chanics at odds with macrorcalism (MAR) [2], a set of 
postulates distilling our common assumptions about the 
macroscopic world. The paradoxes include the violation 
of Bell inequalities in time, also known as Leggett-Garg 
inequalities (LGIs) [2], the non-classicality of weak val- 
ues [3] , and the three-box problem [4] . 

Steady developments in control of single systems for 
quantum information processing have opened the road 
to fundamental investigations of back-action with pho- 
tons [5-7], superconducting circuits [8], and semicon- 
ductor spins [9-11]. Experiments decidedly favor quan- 
tum mechanics over MAR, although loopholes exist in 
each case. Moving beyond fundamental investigation, the 
emergent field of quantum feedback control [12] balances 
the tradeoff between information gain and back-action, 
finding applications in quantum error correction [13], 
qubit stabilization [14, 15], and state discrimination [16], 
for example, where partial measurements are preferred 
over maximally-disturbing projective ones. 

In superconducting circuits, variable-strength mea- 
surement was first demonstrated for a Josephson phase 
qubit [17]. Although destructive for the qubit for one of 
two measurement outcomes, the method allowed prob- 
abilistic wavefunction uncoUapse [18] by two sequential 
partial measurements, firmly demonstrating that back- 
action is phase-coherent [19]. Very recently [20], partial 



measurement of a transmon qubit was realized in cir- 
cuit quantum electrodynamics (cQED) [21, 22] by prob- 
ing transmission through a dispersively-coupled cavity. 
In this case the measurement strength was controlled 
through the number of photons in the probe signal. 

In this Letter, we demonstrate a non-destructive, 
variable-strength measurement of a transmon qubit that 
is based on controlled interaction with an ancilla qubit 
and projective ancilla measurement [Fig. 1(b)]. The key 
advantage of ancilla-based indirect measurement is the 
possibility to accurately tailor the measurement by con- 
trol of the interaction step and choice of ancilla mea- 
surement basis. The kick-back of variable-strength mea- 
surements on the qubit is investigated by conditioning 
tomographic qubit measurements on the result of ancilla 
measurements in different bases, showing close agreement 
with theory. By combining partial and projective mea- 
surements, we observe non-classical weak values, LGI vi- 
olations, and show their predicted correspondence [7, 23] 
in a solid-state system. The ancilla-based indirect mea- 
surement here demonstrated will be extendable to re- 
alization of four-qubit parity measurements needed for 
modern error correction [24]. 

Our planar cQED device [26] consists of two trans- 
mon qubits {Qi and Q2) coupling jointly to a bus res- 
onator (B) and separately to dedicated resonators {Hi 
and H2) used for projective readout [Fig. l(a)]. Lo- 
cal flux-bias lines allow control of individual qubit tran- 
sition frequencies on nanosecond timescales. All mi- 
crowave pulses for individual qubit control and read- 
out are applied through a common feedline coupling 
to Hi and H2. Qubit Qi (Q2) has charging energy 
E'ci{2) — 300 (330) MHz and maximum transition fre- 
quency wSrf2)/27r = 7.37 (6.55) GHz. The bus B has 
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FIG. 1. (color online), (a) Two-transmon, three-resonator 
cQED processor. Quarter-wave resonators Hi and H2 allow 
individual readout of qubits Q\ and Q2 via a common feed- 
line. A resonator bus B (single-photon quality factor 210, 000) 
couples to both qubits. Local flux-bias lines (ports 3 and 4) 
allow independent tuning of qubit transition frequencies with 
~ 1 GHz bandwidth [25]. (b) Scheme for two-step indirect 
measurement of one qubit (Q) through partial entanglement 
with an ancilla qubit [A) followed by projective measurement 
of A. (c) Bloch-sphere illustration of the evolution of A during 
the interaction step, for Q in |0q) and |1q). 



fundamental frequency ujb/'2-'t^ = 6.15 GHz and intrin- 
sic linewidth kb/'^.'k = 30 kHz, and couples to both 
qubits with gB/'i.i: = 36 MHz. Resonant swaps be- 
tween either qubit and B are realized in Tg = 7 ns. 
The readout resonators have fundamental frequencies 
Wj:/i(2)/27r = 7.13 (7.24) GHz and coupling-limited 
linewidth K^f/27r = 9 MHz, and couple to their respec- 
tive qubit with gn/'^'^ = 92 MHz. Projective dispersive 
readout of Qi (Q2) is performed by pulsed homodyne 
detection of feedline transmission at the resonance fre- 
quency of Hi {H2) with Qi {Q2) in the ground state. 
Digitization at the optimal threshold gives 85% (94%) 
single-shot readout fidelity, respectively. 

The interaction step in the indirect measurement is a 
y rotation of the ancilla {A = Q2) by ±e/2, with positive 
(negative) sign for Q Qi in |0q) (|1q)) [Fig. 1(c)]. 
The angle e sets the measurement strength. Note that 
e — 180° makes the measurement projective, as in this 
case A evolves to orthogonal states for |0q) and |1q). 
The (^"dependent y rotation of A is achieved by dressing 
a controUed-z rotation with pre- and post-rotations on 
A [Fig. 2(a)]. The controlled-^ rotation is a three-step 
process: a resonant swap transferring the state of Q to 
B, a photon-controlled z rotation of A, and a resonant 
swap from B back to Q. The acquired two-qubit phase 
SIxaI^w = e is calibrated by varying the wait time Tw and 
the Stark shift on A induced by one photon in B [26]. 
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FIG. 2. (color online). (a) Pulse sequence realizing 
and testing the indirect measurement [Fig. 1(b)] of qubit 
Q = Qi with ancilla A = Q2- Q is first prepared in 
state |6l) = cos(e'/2) |0q) -f sin(6l/2) |1q). (b-d) Ensemble- 
averaged ancilla measurement in the (b) Xa, (c) Ya, and 

(d) Za bases, achieved using pre-measurement rotation 
Ui = Ryj^{—iT /2), identity, and Rx^{~'k/2), respectively. 

(e) Partial-measurement back-action reduces the contrast in 
the final projective measurement of Q. Ideally, (^q) ~ 
(Zq) cos{e/2), independent of Ma basis {Xa here, see [26] 
for the three bases). Inset: Parametric plot of oscillation am- 
plitudes AX'q and AX'a. AXq (AX^) decreases (increases) 
as e increases in the range [0°,180°]. Solid (dashed) curves 
correspond to the model with (without) decoherence during 
the gate. 



Single-qubit phases are nulled by phase-shifting all qubit 
microwave drives after flux pulsing, realizing virtual z- 
gates [27, 28]. 

We characterize the interaction step by performing 
tomographic measurements of Q and A after the in- 
teraction, with Q nominally prepared in the superpo- 
sition state IOq) = cos(6l/2) |0q) sin(6'/2) |1q) and A 
and B in the ground state. As in this part we focus 
purely on the interaction, we correct for readout errors 
using standard calibration procedures [26, 29]. Ideally, 
{X'^) = (ZQ)sin(e/2), {Y^} = and {Z'^} = cos(e/2) 
[(un-)primed notation denotes the (pre-) post-interaction 
state]. These dependencies are well reproduced in the 
data for all choices of e [Figs. 2(b-d)]. Measuring in 
either the Ya or Za basis yields no information about 
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the initial state of Q, as expected. We also measure the 
post-interaction contrast AXg of the Q Rabi oscillation 
[Fig. 2(e)] and compare it to the contrast AX^ of the 
interaction-induced oscillation in A. As e increases, AXg 
increases. Ideally, the quadrature 

(AX'^y + (AX;^)2 = 1 for any e. We observe 

a monotonic decrease from 0.82 at e = 12° to 0.72 at 
e = 180° [Fig. 2(e) inset]. This decrease is reproduced 
by a master equation simulation that includes 4% resid- 
ual excitation and measured decoherence rates for each 
element [26]. 

We now investigate the quantum kick-back of par- 
tial measurements by performing partial state tomog- 
raphy of Q conditioned on the result of the ancilla 
measurement Ma = ±1 in different bases. Results in 
Fig. 3 show the partial-measurement induced distortion 
of a Rabi oscillation of Q for e = 45° and 90° (see 
SOM [26] for other e values) and measurement of Xa, 
Ya, and Za- For Xa, Ma = ±1 kicks Q toward the 
north (south) pole of the Bloch sphere. Ideally, the 
Bloch vector polar angle transforms as 6 — > 9', with 
tan(6l72) = tan(6'/2)[l-M^ tan(e/4)]/[H-AfA tan(e/4)], 
while the azimuthal angle is conserved. Readout errors 
decrease the amplitude of the conditioned curves. The 
difference in the amplitudes is due to asymmetric read- 
out errors [26], taken into account in the model. When 
measuring Ya, conditioning does not distort the Rabi os- 
cillation. This is because the kick-back of Ma = ±1 is a 2; 
rotation of Q by ±e/2, leading to the same x projection. 
Conditioning on a Za measurement produces the most 
striking difference: while Ma = +1 imposes no kick- 
back, AIa = — 1 imposes a z rotation of tt. Ideally, both 
curves are unit-amplitude sinusoids with opposite phase, 
independent of e. However, for e = 45°, the Ma = — 1 
set is dominated by false negatives. As e increases, true 
Ma = — 1 counts become more abundant and we observe 
the expected sign reversal in the conditioned curve with 
e = 90°. Note that despite the difference in conditioned 
curves for the different A measurement bases, the three 
unconditioned curves are nearly identical (See [26] for 
more values of e). This is consistent with the expecta- 
tion that measurement induced dephasing is independent 
of the ancilla measurement basis [12]. 

As a benchmark of the complete indirect-measurement 
scheme, we extract quantum efficiencies rji characterizing 
the loss of quantum information [30] for measurement 
outcome Ma = i- 
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Loss originates in the single-shot readout infidelity of A, 
the residual excitation in A and B, and decoherence of Q, 
A, and B during the interaction. Without decoherence, 
rji would be independent of input qubit state (see [26] 
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FIG. 3. (color online). Partial-measurement kick-back. The 
kick-back on qubit Q induced by partial measurement de- 
pends on the interaction strength e, the ancilla measurement 
basis, and measurement result Ma- Left and right panels cor- 
respond to e = 45° and 90°, respectively. (a,b) Conditioning 
on the result Ma of ancilla measurement in the Xa basis re- 
veals distorted Rabi oscillations of Q. A positive (negative) 
result retards (advances) the oscillation for 9 £ [0°, 180°] and 
advances (retards) it for 8 G [0°,360°]. (c,d) Distortions for 
measurement in the Ya basis. In this case, the kick-back is 
a z rotation by ±e/2, causing an identical reduction of con- 
trast in the conditioned Rabi oscillations. (e,f) Distortions 
for measurement in the Za basis. Ideally, Ma = -f 1 has no 
kick-back on Q, while Ma = — 1 causes a z rotation of tt. The 
lower contrast observed for Ma = — 1 is due to readout errors 
(see text and [26] for details). 



for full expressions) . Using the calibrated infidelities and 
residual excitations, we estimate ri±i = 0.94 (0.94) at 
e = 45° and 0.85 (0.83) at 90°. Including decoherence 
and averaging over the qubit Bloch sphere, we estimate 
7?±i = 0.77 (0.71) at e = 45° and 0.69 (0.60) at 90°. 

Finally, we combine the abilities to perform partial and 
projective measurements to observe non-classical weak 
values, detect LGI violations, and demonstrate their 
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FIG. 4. (color online). Observation of non-classical weak 
values and Leggett-Garg inequality violations (measurement 
strength e = 45°). (a) Pulse sequence. Note that qubit roles 
are swapped (Q = Qi and A = Q\) compared to previous 
figures in order to minimize errors when conditioning on Mq . 
(b) Polar-angle dependence of modified weak value Wm (see 
text for definition and normalization procedure). The MAR 
bound \ Wm\ < 1 is amply exceeded. From the quantum per- 
spective, the extrema in Wm at I and II can be understood 
using the Bloch spheres in (c) and (d), respectively. Ideally, 
at I (II), the Ma = — 1(-|-1) kick-back aligns Q with the +xq 
axis, perfectly correlating Mq — —1 with Ma ~ -|-1(— 1). 
(e) Measured averaged Leggett-Garg operators B± defined in 
text. One of the inequalities — 3 < B± < 1 is violated when- 
ever non-classical Wm is observed. 



correspondence [7, 23] (Fig. 4). The partial measure- 
ment of Q — Q2 is performed via A — Qi and Hi 
and the projective measurement via i?2- We measure 
the partial-measurement average conditioned on strong- 
measurement result Mq = — 1, 



{Ma)\Mq=~i, 



where Ma = {Ma — tIos) /"^pk is the offset and rescaled 
partial-measurement result so that (Ma) = ±1 for 
|0q) (|1q)). mar constrains |Wni| < 1, while quantum 
mechanics allows iV^ml < l/sin(e/2). We call Wm amod- 
ified weak value because it differs in the ideal quantum 
setting (perfect interaction and measurements) from the 
standard definition [3] of the weak value W of operator 



Zq between initial state |0q) and final state |— 7r/2Q), 
-V2q|Zq|0q) 



W : 



Specifically, the digital character of ancilla-based mea- 
surement [26, 31] regularizes Wm near 9 = 7r/2, where W 
diverges [Fig. 4(b)]. In parallel, we consider the averaged 
Leggett-Garg operators 



^4 



±{Ma) T {MaMq) + {Mq) 



Under MAR, a weak, non-invasive measurement and sub- 
sequent strong measurement of a two-level system satisfy 
the LGIs —3 < B± < 1. For the sequence in Fig. 4(a), 
quantum mechanics allows |_B±| < (cos(e) -I- 3) /2. 

Wc capture the smooth crossing of MAR bounds for 
Wm and B± by performing the experiment in Fig. 4(a) 
over a range of initial qubit states \0q). We observe 
a maximum Wm of 1.57 ± 0.08. In turn, the averaged 
Leggett-Garg operators and B- peak at 1.25 ± 0.05 
and 1.19 ±0.06, respectively. The data clearly show that 
one of the two LGIs is violated whenever Wm is non- 
classical. This correspondence, predicted in Ref. 23 and 
previously only demonstrated with photons [7] , becomes 
the more interesting upon noting that B± averages all 
measurements while Wm uses only the post-selected frac- 
tion for which Mq = — 1. 

In conclusion, we have realized an indirect measure- 
ment of a transmon qubit with high quantum efficiency 
and tunable measurement strength. Our scheme consists 
of a partially entangling interaction between the qubit 
and an ancilla, followed by projective ancilla measure- 
ment using a dedicated, dispersively-coupled resonator. 
We have measured the kick-back of such measurements 
on the qubit as a function of interaction strength and 
ancilla measurement basis, finding close agreement with 
theory. Non-classical weak values are observed upon 
conditioning ancilla measurements on the outcome of a 
projective measurement of the qubit. Their predicted 
correspondence with LGI violations is demonstrated for 
the first time in a solid-state system. The combination 
of high-quality factor bus, individual readout resonators 
and feedline here demonstrated constitutes a scalable ar- 
chitecture [32] with frequency-multiplexable single-qubit 
control and readout [33]. Future experiments using this 
architecture will target the realization of an ancilla-based 
4-qubit parity measurement as needed for surface-code 
quantum error correction [24]. 
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DEVICE FABRICATION 

The two-qubit, three- resonator chip was fabricated on 
a sapphire substrate (430 fim thick, C-Plane). Follow- 
ing in-situ cleaning of the substrate in Ar for 2 min, a 
65 nni thick niobium titanium nitride (NbTiN) film [1] 
was DC-sputtered. Superconducting coplanar waveguide 
structures (10 /xm central conductor width, 4.2 /um gaps) 
were then defined using a negative electron-beam resist 
(SAL-601) and reactive-ion etching in a SF6/O2 plasma. 
Finally, the two transmons were patterned by standard 
electron-beam lithography and Al double-angle evapora- 
tion (20 nm bottom and 70 nm top layer thicknesses, with 
intermediate oxidation for 10 min at 0.55 mBar). During 
fabrication, the device was exposed three times to an O2 
plasma to remove organic residues: before evaporation, 
after lift-off, and after dicing the sample to 2 mm x 7 mm. 

EXPERIMENTAL METHODS 

This section consists of four fuUy-captioned figures pro- 
viding further detail on experimental methods. The com- 
plete setup, both inside and outside the dilution refriger- 
ator, is illustrated in Fig. SI. The detailed flux-pulsing 
scheme showing the bias points chosen for single-qubit 
control, readout and interactions with the bus is shown 
in Fig. S2. Figure S3 presents the calibration of the in- 
teraction step for the values of measurement strength e 
used. Finally, Fig. S4 demonstrates the use of resonators 
Hi and H2 for individual readout of Qi and Q2, respec- 
tively, and the extraction of single-shot readout fidelities 
and residual qubit excitations. 

EXTENDED RESULTS 

This section consists of four fully-captioned figures 

extending the results and backing claims in the main 
text. Figure S5 demonstrates that partial-measurement 
induced dephasing is independent of the basis chosen 
for the ancilla measurement. Figure S6 shows the 
measurement-induced back-action at more values of mea- 



surement strength e than Fig. 3. Figure S7 shows the 
raw data contributing to the measurement of weak val- 
ues, complementing Fig. 4(b) and showing raw measure- 
ments at more values of e. Figure S8 shows the three 
terms contributing to the averaged Leggett-Garg opera- 
tors B±, complementing Fig. 4(c). 



THEORY 
Hamiltonian model 

We consider a system of two transmons (Qi and Q2) 
coupled to one bus resonator. The system Hamiltonian 
is 

2 

i=l 

where Hb and Hqi describe the non-interacting dynam- 
ics of the bus and qubit Qi, respectively [2]: 

Hb = nujBO'a, Hq^ = ^huij^Qi \jQi) {jQi\ . (S2) 

j 

Here, ujb is the bus resonance frequency and fkOj^Qi is 
the energy of the j^^ level of qubit Qi. In the transmon 
regime Ej/Ec ^ 1 valid here [2], the interaction between 
a transmon and the bus can be modeled by an extended 
Jaynes-Cummings-type coupling [2]: 

^Qi = fi^9j+i,j {\j + Igi) a + h.c.) , (S3) 
j 

with coupling strengths gj+ij = \J j + IffO) where 2gQ is 
the vacuum Rabi splitting. The driving Hamiltonian is 

^Qi(i) =E.,,,d (Q(t)x/J+T|i + lQi)0-Qi|e-'-^*+^<* 
+ h.c.) , 

(S4) 

where ed(t), <f>d,, and denote the amplitude, phase, and 
frequency of drive d, respectively. When each qubit is far 
detuned from the bus, i.e. (?j+ij/|u;j+i,Qi— Wj,Qj— wb| <C 
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1, the system is described by a dispersive Hamiltonian [2] 
with transmon and transmon-bus interaction terms 

^Qi = ^^^3,Qi \jQi) ijQil + Xj,Qi \j + lQi> {j + iQil , 



= -fiXo,Qia^a\OQi) (OQj| 

+ fi'^{Xj-i,Qi - Xj,Qi)a^a \jQi) {jQil ■ 



(S5) 



Here, Xj,Qi = fff+i - Wj.gi - ws) is the dis- 
persive bus-transmon couphng. The dispersive regime 
allows for an analytical solution for the shift of all res- 
onance frequencies, and simplifies the modeling of the 
gate sequence. We therefore use Eq. (S5) to simulate the 
dynamics in the off-resonant passages with Hgi — )• Hq^ 
and Iq, ^ in Eq. (SI). 

Apart from the coherent dynamics, the coupling of the 
system to the environment leads to dissipative evolution. 
Assuming weak coupling, the Markovian master equation 
describing the system evolution is [3] 



p = -i[H, p] + kD [a] p -h k'V [a)] p 

2 

1=1 

2 

i=l 

+ Y.^-^V[Zq,]p^Cp, 



(86) 



Partial measurement : ideal CEise 

For reference, we analyze the partial-measurement 

scheme with perfect interaction and measurement steps. 
In the interaction step, the qubit and ancilla evolve as 



with 



U = cos{e/4)lQlA - ism{e/A)ZQYA. 



Upon performing a projective measurement of the ancilla 
with operator Oa = J2i N^) (unnormalized) 
post-measurement qubit state for result Ma = A, is 

|*'Q) = fi,J*g), 

where fl\^ = {ia] U \0a)- The probability of getting mea- 
surement result Ma = Aj is 

Pm.=a. = (*q| 4.^a. |*q) • 

When the measurement results are disregarded, the post- 
measurement qubit density matrix is 

i 

= TrA [C/Pq|Oa) (Oa|C/^] 

= cos2(e/4)pQ + sin2(e/4)ZQPQZQ, 



where V[L\p= {2LpL^ - L'^Lp - pLL^)/2, k {k') is the 
bus photon damping (excitation) rate, and 7oi, Qi, 7io, Qi) 
and 70, Qi are relaxation, excitation, and pure dcphas- 
ing rates for Qj. Here, we have explicitly truncated the 
Hilbert space of each transmon to the lowest two lev- 
els. Equation (S6) is conveniently solved in Liouville 
space, where p is a vector. The solution is given by 
p{t) = e^*/9(0), which can be solved by numerically diago- 
nalizing the propagator C. In the simulation, we truncate 
the bus Hilbert space to n = {0, 1, 2} photons. 

The model is fit to the data with the following fit pa- 
rameters: the frequency of Qi during the interaction with 
B, the two single-qubit phases, and the amplitude error 
in the pre- measurement rotation for Q. These parame- 
ters are independently fit for each pair of simultaneous 
measurements of Qi and Q2- For conditioned measure- 
ments, we also fit the readout fidelities of the conditioning 
readout. Fixed model parameters include 4% excitation 
in each element [measured for Q2 (Fig. S4) and assumed 
equal for Qi and B], the measured energy-equilibration 
times 1.4, 2.5, and 5.3 ps for Qi, Q2, and B, respectively, 
and pure-dephasing time 1.6 (1.8) ps for Qi {Q2)- Typi- 
cal best-fit values give an absolute error within 5° for e, 
10°(3°) for the phase of Qi {Q2), and 2% for the rotation 
amplitude. 



where pq = I^'q) (^^qI is the initial qubit density matrix. 
Clearly, pg is independent of the choice of ancilla mea- 
surement basis. It is straightforward to show that the 
transformation 




of the qubit Bloch vector is 




cos(e/2) (Xq) 
cos(e/2) (Yq) 
{Zq) 



Ancilla measurement in Xa 



The operation elements evaluate to 

^±1 = ^ cos(e/4)/Q ± sin(e/4)ZQ, 
giving the measurement result probabilities 

^M.=±i = ^(l±sin(e/2)(ZQ)) 
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and the measurement average 



and thus a measurement average 



{Ma) = {Zq) sm{e/2). 

We briefly visuaHze how these operation elements kick 

a qubit initially in the pure state \iI>q) — \0t4>) = 
cos(6'/2) |0q) + e''^sin(6'/2) Because these opera- 

tion elements are real-valued and diagonal, the post- 
measurement state |V'q) — Wi4'') has the same az- 
imuthal angle, </>' = (f). The polar angle transforms as 



tan(6'72) =tan(6'/2) 



cos(e/4) - MAsin(e/4) 
cos(e/4) + MAsin(e/4)' 



We consider a few special cases. For a qubit initially 
on the equator of the Bloch sphere {9 = it/ 2), a posi- 
tive (negative) measurement kicks the qubit toward the 
north (south) pole, decreasing (increasing) its polar angle 
by e/2. A qubit initially at one of the poles (^ = 0, tt) re- 
mains at the pole regardless of measurement result. For 
all other cases, the update formula shows that the change 
in polar angle is not equal in magnitude. In the north- 
ern (southern) hemisphere, a positive result decreases the 
polar angle less (more) than a negative result increases 
it. 



Measurement in Ya 



The operation elements are 

fl±i = ^(cos(e/4)T«sin(e/4)ZQ) 
= 7I^^Q(±e/2), 



yielding 



2' 



(Ma) = 0. 



The kick-back of Ma = ±1 is evidently a z-axis rotation 
of Q by ±e/2. This kick-back is independent of the initial 
qubit state. 



Measurement in Za 



(Ma) = cos(e/2). 

For the most likely result {Ma = +1), there is no back- 
action. The rare result {Ma = — 1) rotates the qubit by 
TT around the z axis. 

Modified weak value: ideal CEise 

Consider performing the general partial plus projective 
measurement scheme (Fig. 4) starting with the qubit in 
IV'q) = 1^) and the ancilla in \0a)- A projective measure- 
ment of Xq on Q following the interaction step has for 
operation elements on A 

OM«=±i = (±7r/2Q|[/|^Q)|0^). 

The (unnormalized) post-measurement ancilla states are 

I^'a) =nM^=±i \0a) 

= cos{e/4){±n/2Q\9Q)\0A) 

+ sin(e/4) (±7r/2Q|ZQ|^Q)|U). 

From these we can calculate the expectation value of Ma 
{Xa basis) conditioned on Mq = ±1: 

/M \ — sin(e/2) 

( a)\Mq=±i -cos2(e/4)+sin2(e/4)|M^(ZQ,±7r/2,^)|2^ 



Re{W{ZQ,±n/2,9)}, 



where 



W{Aq,9',9) = 



{0q\Aq\9q) 
{0'q\Oq) 



is the weak value of Hermitian operator Aq between ini- 
tial state \9q) and final state \9'q). 



The operation elements are 



yielding 



n+i = cos(e/4)/Q 
n_i =sin(e/4)ZQ, 



Pma=+i = cos2(e/4) 



Ma 



sin2(e/4). 



Me2isurement model 

In order to include the readout errors in the model 
curves, we numerically calculate the density matrix p of 
the two-qubit, bus-resonator system following the mea- 
surement pre-rotations. Next, we use the calibrated read- 
out errors for Q and A to calculate conditioned and un- 
conditioned averages. For example. 
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IM \ — Pao[C2Fqo-1)pao,qo-(2Fqi-1)pao,qi] + ('^-Fai)[(2Fqo-1)pai,qo-(2Fqi-1)pai,qi] 

V^^QI\Ma = + 1 ~ PAoFAO+PAlil-FAl) ■ 

if./r \ _ F^i[(2FQo-l)p^i,Qo^(2fQi-l)pAi,Qi] + (l-FAo)[(2FQo-l)pAo,Qo-(2FQi-l)p^o,Qi] 

V^-^QI\Ma = -1 ~ PaiFai+pao{1-Fao) ■ 

{Mq) = {2Fqo - 1)PQ0 - {2Fqi - 1)PQU 



where 

PAi = Tr[|zA) {iA\p] , 

PQj = Triljg) (jglp] , 

PAi,Qj = Tv[\jQiA) {jqiAlp] 

are the probabihties of A being in 1?^), Q being in Ijq), 
and A and Q being in {Jqia), respectively, and Fai and 
Fqj are the single-shot readout fidelities for A and Q 
calibrated in Fig. S4. Conditioned and unconditioned 
averages for Ma are similarly calculated. 

Quantum efficiency 

We consider the evolution of the qubit density matrix 
Pq conditioned on specific ancilla-measurement results: 

„ _ f Poo Poi \ MA=i I _ ( Poo P'oi \ 
""-{pro Pn) ''^-[p[o P'nj- 



We define the quantum efficiency rji for measurement out- 
come Ma = i by 

I Poll ^ I Poll 
VPooPu ' \/PooPii ' 

capturing the loss of information about the qubit as a re- 
sult of the partial measurement [4] . For an ideal indirect 
measurement with perfect interaction and readout steps, 
it is straightforward to show rj^i = rj-i = 1. Including 
infidelity in the ancilla readout gives, for Ma = +1, 

P'q OC FaqO+iPqO]^! -hFAiO.iPgfiLi, 

where Fai = 1 — F^i. The expression for Ma = — 1 is 
obtained by substituting Fao — >■ Fao and Fai — >■ -Fai- 
The outcome-specific quantum efficiencies for our choice 
of qubit-ancilla interaction and measurement in the Xa 
basis become 



J 



(FAo+FAi)cos(e/2) 
^J[Fao + Fai + {Fao - Fai) sin(e/2)] [Fao + Fai - {Fao - Fai) sin(e/2)] 

and 

{FAo + FAi)cos{e/2) 

V-i = I _ _ 

^ [Fao + Fai - [Fao - Fai) sin(e/2)] [Fao + Fai + {Fao - Fai) sin(e/2)] 

These quantum efficiencies do not depend on the initial state of the qubit. The asymmetry in the ancilla readout 
fidelities causes ?7+i and ?7_i to differ. Including the residual excitation of the ancilla, PgA, further reduces the 
quantum efficiencies to 

{Fao + Fai) cos{e/2) 

^[Fao + Fai + (1 - 2PeA)(-FAo - Fai) sin(e/2)] [Fao + Fai - (1 - 2PeA)(-FAo - ^ai) sin(e/2)] 

and 

{Fao + Fai) cos{e/2) 

V-i = I _ 

^[Fao + Fai - (1 - 2PeA)(PAo - Pai) sin(e/2)] [P^o + Pai + (1 - 2P^a){Fao - Fai) sin(e/2)] 



Including decoherence of the ancilla, qubit and bus In this case we can average rjj over the surface of the 
makes jj+i and dependent on the qubit input state. qubit Bloch sphere to arrive at single numbers. Lack- 
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ing closed-form formulas, we rely on the master equation 
simulation to calculate 



Vi 



J r]iSm{9)d9d(j) 
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FIG. SI. Experimental setup. Arbitrary waveform generators Tektronix AWG520 and AWG5014, with 10- and 14-bit resolution, 
respectively, and 1 ns sampling rate produce voltages directly applied to the flux-bias lines, the single-sideband I-Q modulation 
envelopes for the microwave tones driving single-qubit rotations, and the pulse envelopes for measurement tones. Flux pulses 
are conditioned by a series combination of DC block, attenuation, LC low-pass fllter and homemade coaxial eccosorb filter 
(based on Ref. 5) before reaching ports 3 and 4. All measurement and qubit-drive pulses are combined at room temperature. 
Inside the dilution refrigerator, they are coupled to the feedline input (port 1) following 50 dB attenuation. On the feedline 
output (port 4), an amplification chain with ~ 100 dB gain, two I-Q mixers and a two-channel averaging digitizer (1 ns, 8-bit 
sampling) process the two readout signals. 
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FIG. S2. Detailed sequence of microwave and flux pulses realizing the indirect measurement in Fig. 2(a). The qubits are 
detuned from each other to minimize crosstalk in qubit control and readout. Qi — Q {Q2 = A) is pulsed to the bias point 
[6.69 (5.74) GHz] where all single-qubit operations are performed. All single-qubit gates are realized with resonant DRAG [6] 
pulses, with standard deviation a — 4 ns, and ±2(7 truncation. The rotation axis is set using I-Q (vector) modulation (see 
Fig. SI). The controUed-z rotation is realized by coherently swapping the Qi state into the bus B in Ts = 7 ns, waiting 
a calibrated time Tw (Fig. S3), and then swapping the B state back onto Qi. The photon-number dependent shift of Q2 
during Tw produces the two-qubit phase that can partially entangle Qi and Q2- The measurement pre-rotations Ui for ancilla 
measurement in the Xa, Ya, and Za bases are J?y^(— 7r/2), identity, and Rx^{—n/2), respectively. Before measurement, Qi 
is flux-pulsed to 6.39 GHz and a tt pulse on the 1-2 transition of Q2 is applied to maximize fidelity. A 2 ns buffer is inserted 
between adjacent pulses to avoid any overlap. 




Calibrated (t>+E/2 (deg) 

FIG. S3. Calibration of the partial-measurement interaction step, (a) Pulse scheme for measuring the single-qubit phase 
of Q2 and the two-qubit phase acquired during fiux pulsing. We measure the shift in the azimuthal phase of the Q2 Bloch 
vector with Qi prepared nominally in |0) (6» = 0) and in |1) (6* = 180°). The peak of the 6* = (180°) curve in (b) [(c)] is 
matched to azimuthal angle tj) = —180° =F e/2. In this way we implement virtual 2-gate corrections [7] as common in nuclear- 
magnetic-resonance experiments [8]. (d) Table of calibrated e for various conditions of the Q2 fiux pulse in between the Qi-B 
swaps. The duration Tw (1 ns resolution) and the flux-pulse amplitude AV2 are used to coarse- and fine-tune e, respectively. 
The single-qubit phase acquired by Qi is calibrated by interchanging the rotations applied to Qi and Q2, and is similarly 
compensated. 
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FIG. S4. Characterization of individual qubit readouts. (a,b) Averaged pulsed measurement of feedline transmission as a 
function of frequency near the fundamental-mode (quarter-wave) frequencies of (a) H\ and (b) H2 immediately following 
preparation of with tt pulses G {0, 1} denote the state of Q2 and Qi, respectively). Here, we use 280 ns integration 
time and —121 dBm (—111 dBm) input power at port 1. Evidently, Hi {H2) is predominantly sensitive to the state of Qi 
(Q2), with -6.4 MHz (-2.0 MHz) dispersive shift. (c,d) Single-shot histograms at fn^ = 7.130 GHz and /^^ = 7.243 GHz 
[600 ns integration time, —110 dBm and —111 dBm input power, respectively]. Using the best fit of a double Gaussian to the 
H2 histograms for Q2 nominally in |0), we estimate that 4% of counts fall in the peak corresponding to |1). We attribute this 
fraction to residual excitation of Q2- This hypothesis is supported by other measurements (not shown) at variable power and 
duration, giving similar fit results. The power dependence observed in histograms for Hi with Qi nominally in |0), however, 
does not allow this analysis. When modeling, we thus assume a residual excitation of Qi equal to that measured for Q2. (e,f) 
Cumulative probability of histograms in (c,d). The readout contrast is I — cq — ci = 0.779±0.005 for Hi and 0.870±0.007% for 
H2. (g,h) Readout error model for Qi and Q2- Accounting for the contrast reduction induced by residual qubit excitation and 
assuming perfect pulses, we extract single-shot readout fidelities -Fb.Qi = 0.93, Ji,qi ~ 0.92, Fo,q2 = 0.99, and F1.Q2 — 0.95. 
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Rabi-rotation angle 6 (deg) 

FIG. S5. {X'q) as a function of the qubit Rabi-rotation angle 9, for three choices of basis for the ancilla measurement Ma- (a) 
Xa, (b) Ya, (c) Za- Panel (a) replicates the data in Fig. 2(e). Clearly, unconditioned measurements of Q do not depend on 
the choice of ancilla measurement basis. 




FIG. S6. Measurement back-action over a wide range of measurement strength e. Panels (b), (c), (g), (h), (1), and (m) replicate 
the data in Fig. 3. 
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FIG. S7. Weak-value measurements over a wide range of measurement strength e. Panel (b) shows the raw data used in 
Fig. 4(b). 
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FIG. S8. The three terms (Ma), (MaMq), and (Mq) contributing to the averaged Leggett-Garg operators B± shown in Fig. 4. 
As discussed in the main text, Ma is the partial-measurement result, offset and rescaled so that (Ma) ~ ±1 for |0q) (|1q))). 
In the ideal quantum setting, the two-qubit correlation term (MaMq) = for all d and the reduced contrast in (Mq) arising 
from partial-measurement kick-back is (Mq) — sin(^) cos(e/2). The vertical asymmetry in (Mq) observed in data and model 
is due to asymmetric errors in the readout of Q (Fig. S4). 



